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Abstract 

We have measured the magnetic field dependence of the ordered antiferromagnetic moment and 
the magnetic excitations in the heavy-fermion superconductor URu2Si2 for fields up to 17 Tesla 
applied along the tetragonal c axis, using neutron scattering. The decrease of the magnetic intensity 
of the tiny moment with increasing field does not follow a simple power law, but shows a clear 
inflection point, indicating that the moment disappears first at the metamagnetic transition at ~40 
T. This suggests that the moment m is connected to a hidden order parameter ip which belongs 
to the same irreducible representation breaking time-reversal symmetry. The magnetic excitation 
gap at the antiferromagnetic zone center Q=(1,0,0) increases continuously with increasing field, 
while that at Q=(1.4,0,0) is nearly constant. This field dependence is opposite to that of the gap 
extracted from specific-heat data. 

PACS numbers: 75.25.+Z, 75.30.Kz, 75.30.Cr, 75.50.Ee, 74.70.Tx 
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One of the most enigmatic problems in heavy-fermion physics is to reconcile the smallness 
of the ordered antiferromagnetic (AFM) moment of URu2Si2 with the huge A-anomaly in 
the specific heat, Cp, at Tjv = 17.5 K. Neutron [0, 0, ||, |, |^ and resonant X-ray [|| scat- 
tering experiments have shown that a weak static antiferromagnetic dipolar order with the 
moment along the c axis of the body-centered tetragonal structure (space group I4/mmm) of 
the heavy-fermion superconductor (Tc = 1.2 K) URu2Si2 is established below T/v. The mag- 
nitude of the moment (0.03 yU^/U-atom) is too small to reconcile with the jump in Cp. The 
specific-heat anomaly suggests the opening of an energy gap over part of the Fermi surface 

|l^. Inelastic neutron-scattering measurements below show well-defined dispersive 
magnetic excitations polarized along the c axis, which are consistent with a singlet- singlet 
crystal-field model |Il|. However, such a model predicts an ordered moment much stronger 



than that observed. Recent reports |Tl|, on that URu2Si2 samples are inhomogeneous 
with a larger moment in small regions are not yet experimentally confirmed and are still a 



matter of controversy and debate |13]. Also, they do not resolve the discrepancy between 
the small moment and the specific-heat anomaly. 

In this Letter we present the magnetic field dependence of the AFM component of the 
ordered moment and of the gaps in the magnetic excitation spectrum for fields H up to 17 
T applied along the c axis. The high magnetic fields now available for neutron scattering 
experiments at the HMI made it possible to observe that the magnetic moment in URu2Si2 
does not vanish at ~15 Tesla as was suggested from earlier measurements using a simple 



power-law extrapolation 0, |I4| . Instead, we observe a clear inflection point of the magnetic 
intensity near a field of 7 T and a finite moment at 17 T. Such a behavior can be obtained 
if there is a linear coupling between a hidden primary order parameter and the magnetic 



moment |jT5[ . Inflection points have also been suggested in other theoretical schemes, which 
will be discussed. The gap A in the magnetic excitation spectrum at the antiferromagnetic 
zone center Q=(1,0,0) increases strongly with H. The second gap V at Q=(1.4,0,0) is nearly 
field independent. 

Magnetic fields up to 17 T at low temperatures were obtained by mounting the sample 
between two dysprosium rods in a 14.5 T superconducting magnet ||16|. A high-quality cylin- 
drical URu2Si2 single crystal of diameter 6 mm and height 3.8 mm was cut by electroerosion 
from a larger well-characterized and annealed crystal already used in previous work WA 
and mounted with the c axis parallel to the vertical field of the magnet. The neutron scat- 
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tering experiment was made on the FLEX triple-axis spectrometer of HMI, using pyrolytic 
graphite as vertically focusing monochromator and horizontally focusing analyzer. Measure- 
ment of the magnetic moment was done using the spectrometer in W (RLR) configuration 
with a wave vector of 1.48 A^^ To remove the higher-order harmonic contribution from the 
monochromator, we used a 10-cm thick liquid-nitrogen cooled Be filter in the incident beam 
and a graphite filter in the scattered beam, the latter oriented to reflect out the 2nd order 
harmonics. The inelastic measurements were done with the spectrometer in a "long-chair" 
(RRL) configuration, a fixed final wave vector of 1.55 A^^, and the Be filter in the scattered 
beam. The energy resolution at elastic energy transfer was 0.14 meV at full width at half 
maximum. Most of the elastic and inelastic data were collected at a temperature of T = 2 
K for magnetic fields H between and 17 T. 

The integrated intensity of the AFM Bragg peak, which is proportional to the square 
of the ordered magnetic moment, m?, was measured by rocking scans at Q=(1,0,0), as 
illustrated in Fig. [l|. Measurements above T/v (at high fields) show clearly the absence 
of any second-order contamination. The finite correlation length is field independent, as 
reported in earlier measurements Figure ^ shows the integrated Bragg peak intensity 



as a function of applied field H \\ c. The figure combines the present data from HMI with 
earlier measurements at fields up to 12 T from the IN14 spectrometer at the ILL [Q, using 
the field points at and 8 T to scale the results. The magnetic moment decreases with 
increasing field, but clearly deviates at higher fields from the standard expression 

=ml\l- [H/Hmy] (1) 

with 7 = 2. Equation (P describes the magnetic field dependence at T = of a (single) 
order parameter m (the magnetic moment) for a second-order phase transition, and is easily 
derived in Ginzburg-Landau theory using a free energy of the form 

F = a{h'^ - t)m^ + /3m^ (2) 

with if: = 1 — T jTfq and h = H/Hm, where T/v and Hm are the critical temperature and field 
for the magnetic order, respectively, and a and (3 are constants. The critical field Hm was 
extrapolated to ~15 T in earlier low-field work |14|. Modifying the 7 exponent (a value 
of 3/2 was used in Ref. 0]) improves the fit but does not remove the disagreement with the 
data, which shows a clear inflection point at ~7 T. 



An inflection point in m{H) was obtained in a model proposed by Shah et al. where 



the magnetic moment m is a secondary order parameter, which is coupled linearly to a 
hidden primary order parameter ip, which breaks time reversal symmetry. The magnetic 
intensity can in this case be written as 

where He is the metamagnetic field above which the heavy-fermion state is suppressed 
{He =35. 8-39. 4 T l]!^; we use Hc=4:0 T) and 5 is a parameter. Shah et al. introduces a 



characteristic field for the magnetic moment through Hm = Hc/{1 + 25)^^"^. Equation 
(I) gives an excellent description of the observed field dependence, as illustrated in Fig. |^, 
with Hm = 10.4 ± 0.2 T as the only adjustable parameter. The linear coupling term oc rmjj 



was originally suggested in Ref. 0. Symmetry arguments |jT5|, show that ip must have 
the same symmetry as m, i.e. ip belongs to the irreducible representation (IR). The 
physical interpretation of is either in terms of triple-spin correlators (octupolar moments 
or couplings between dipolar and/or quadrupolar moments) or a dipolar moment of the 
conduction electrons, which have such a rapidly decreasing form factor that the magnetic 
intensity is negligibly small already at the first observable magnetic Bragg peak [|T^ . 

The present finding of an infiection point in m{H) rules out the possibility that the pri- 
mary order parameter is even under time reversal. Taking the simplest coupling term m^?/'^ 
with the right AFM symmetry |jT5|, |18|, the field dependence of m is the same as in Eq. 



151 , and hence in contradiction with the present results. Recent ^^Si NMR data suggest 
in fact that ip breaks time- reversal symmetry UTS). The existence of an infiection point in 
m{H) also excludes the scenario with two decoupled order parameters and accidentally close 
transition temperatures [jl9|, since m{H) would then follow Eq. (|l]). 



While a hidden order parameter that breaks time-reversal symmetry is the most likely 
scenario, an infiection point in m{H) can also be obtained in a strong coupling model with 



the magnetic moment as the only order parameter If the total moment changes between 
the normal and the magnetically ordered state due to an increase of the magnetic coupling 
below T/v, then the (3 term will depend on the external field h and is replaced by (3{1 + r]h'^) 
in Eq. (0), which leads to 

1 - /?2 

Equation (^) provides a good description of the present neutron data, with rj = 11. A ± 0.7, 



as seen in Fig. ^ However, the model does not resolve the discrepancy between the small 
moment and the specific-heat anomaly. 

Okuno and Miyake |21[] have proposed a model which combines the localized character 
of a the singlet-singlet crystal-field scheme with an itinerant character of the 5/ electrons. 
They assume further that there is a perfect nesting at the AFM zone center Q=(1,0,0). For a 
suitable parameter set, they obtain an infiection point in the field dependence of the magnetic 
moment whose shape resembles that of the present results. However, for the parameter sets 
they employ, the inflection point occurs at magnetic fields several times higher than the 
observed inflection point and it seems difficult to obtain a semi-quantitative agreement with 
experimental values of the size of the ordered moment, the transition temperature, and the 
inflection point, simultaneously. 

We note in passing that a magnetic fleld applied along the c axis induces a localized 
ferromagnetic component, resulting in a ferrimagnetic structure with two different moments 
on the uranium-ion sites. The magnitude of the ferromagnetic moment at the metamagnetic 
transition can be estimated to 0.32 fiB, based on an extrapolation of low-fleld polarized 



neutron diffraction measurements 22 



The magnetic excitations at the local minima Q=(1,0,0) and (1.4,0,0) of the dispersion 
curve have been measured at low temperature (T = 2 K) for magnetic flelds up to 17 T using 
the FLEX spectrometer of HMI. These data have been combined with earlier measurements 



at i7 < 12 T using the IN14 spectrometer of ILL \T^]. The results are shown in Fig. |^. 
The measurements on IN14 have better statistics, due to higher flux and a larger sample 
volume (the sample volume is very limited in the 17 T magnet). Instrumental resolution 
effects in combination with the steep dispersion introduce an asymmetric high-energy tail of 
the sharp excitation at the magnetic zone center Q=(1,0,0). In such a case, the excitation 
energy can be extracted by fltting the convolution of a Gaussian function, which describes the 
energy resolution, with an exponential saw tooth function exp[— (£" — A)/e], which describes 
how the higher-energy excitations are coupled into the measured spectrum via the flnite Q 
resolution. The fleld dependence of the parameter e, which has no physical meaning, was 
determined from measurements at high energy transfers on IN14, and extrapolated to higher 
fleld values. This method gives a very accurate description of the measured spectra, as shown 
in Fig. 0(a), and the obtained excitation energy A lies, as expected, at approximately the 
half-height of the leading (low-energy) edge of the measured intensity. The peak shape of 
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the excitation at the second minima, Q=(1.4,0,0), is symmetric [see Fig. |5(b)] because of the 
smaller dispersion, and a Gaussian fit is sufficient to extract the corresponding excitation 
energy V. 

The magnetic field dependence of the gap energies A and V are shown in Fig. ^ The 
gap at the AFM zone center increases rather strongly with applied field, without any sign 
of saturation at higher fields. It is well described by 

A{H) = [Al + {bHr]^' (5) 

with Ao = 1.59 (2) meV and b = 0.114 (1) meV/T. The gap V = 4.51 (2) meV at 
Q= (1.4,0,0) is nearly constant, with only a slight linear decrease in energy with increasing 
field, with slope 5.4 ± 2.7 /xeV/T. It is interesting to note that the high-field extrapolation 
of the energy gaps A (if) and V{H) crosses at a field close to the metamagnetic transition 
He- This is very interesting in the light of the appearance of a new high-field phase above 



36.1 T p3|, which could be due to Fermi surface effects requiring a different nesting vector, 
such as the one corresponding to Q= (1.4,0,0). 

Specific heat |Ty] and electrical resistivity data taken with a magnetic field if || c up 



to 17.5 and 16 T, respectively, have been analyzed in terms of a field-dependent energy gap 
related to the excitation gap observed in neutron scattering measurements. In both cases, 
the Neel temperature and the gap energy decrease quadratically with the field, following the 
relations r^(ii) = Tn{0)[1 - [H/H^f] and A(ii) = A(0)[1 - {H/Ho)% respectively, with 
iio being close to the metamagnetic field. Very recent specific-heat measurements at fields 
up to 45 T pS) show that T/v vanishes at a critical field of 35.9 T. The field dependence of 
T/v is essentially quadratic. Thermal expansion data with a magnetic field ii || c up to 
25 T have been fitted by the empirical function Aexp{—A/T) for the magnetic contribution. 
A quadratic decrease of A and Tn was found with an extrapolated critical field of 37 T and 
a gap energy of 115 K at zero field. 

The important point is that the energy gap extracted from different bulk measurements 
has the opposite magnetic field dependence to the gap A observed directly by neutron scat- 
tering at the AFM zone center. This low-energy gap actually increases strongly with field. 
The second gap V shows a very slight linear decrease, much slower than a quadratic de- 
crease observed in macroscopic measurements. While specific heat and electrical resistivity 
measure some average gap energy, it is difficult to imagine how an average of a strongly 
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increasing low-energy gap and a nearly field- independent high-energy gap would give rise 
to a decreasing gap. It is therefore clear that the gap from bulk measurements is unrelated 
to the gap(s) in the magnetic excitation spectra. This is in stark contrast to the transition 
temperature, which has the same field dependence in specific heat [ITUl ESI |23|, electrical 



resistivity thermal expansion [^], and neutron scattering ^ measurements. 

Unfortunately, due to the difficulty to measure T/v accurately for small moments with neu- 
tron scattering, there are no systematic neutron measurements of T/v (if). It is interesting 
to note that the height of the anomalies in the specific heat ACp [l^ and resistivity Ap |24 
as well as the intensity of the magnetic gaps A and V (see Fig. H) are field independent. 

In conclusion, our neutron scattering measurements on the heavy-fermion superconduc- 
tor URu2Si2 show that the magnetic field dependence of the antiferromagnetically ordered 
moment m has a clear inflection point at ~ 7 T and remains finite even at 17 T for fields 
applied along the c axis. This strongly suggests the existence of a hidden order param- 
eter ip that breaks time-reversal symmetry. Symmetry arguments show that the hidden 
order parameter belongs to the same irreducible representation ^423 as m and corresponds to 
(ordinary) dipolar moments of the conduction electrons or triple-spin correlators of the 5/ 
electrons. Our neutron scattering measurements of the (gapped) magnetic excitations show 
that they are not related to the gap extracted from e.g. specific heat measurements, since 
the magnetic field dependence of these gaps is completely different. Instead, it is the hidden 
order parameter that gives rise to the large anomaly in specific heat and other macroscopic 
properties. The extrapolated crossing of the two magnetic gaps might explain the new phase 
recently discovered above the metamagnetic transition. 

We acknowledge fruitful discussions with G. Amoretti, N.H. van Dijk, J. Flouquet, I. 
Fomin, V.P. Mineev, K. Miyake, P. Santini, J. Schweizer, M.B. Walker, and M.E. Zhit- 
omirsky. A. A. Menovsky, P. Lejay, and A.D. Huxley contributed to sample growth and 
preparation. We thank M. Jaime for giving us access to his high-field specific-heat data 
prior to publication. 
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FIG. 1: Rocking scans of the AFM Bragg peak for different values of the magnetic field. The lines 
are Gaussian fits. 
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FIG. 2: Integrated intensity of the AFM Bragg peak as a function of magnetic field. Squares 
(circles) are from FLEX (IN14). The dotted, solid, and dashed lines correspond to fits given by 
Eqs. (|l]), (^), and (^), respectively. 
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FIG. 3: Energy scans of the magnetic excitations at the two local minima Q=(1,0,0) (a) and 
Q=(1.4,0,0) (b) for magnetic fields of 0, 12, and 17 T. The intensity scales for the data on the 
IN14 and FLEX spectrometers are scaled to coincide. The lines are fits described in the text. 
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FIG. 4: The gap energies at Q=(1,0,0) (solid symbols) and Q=(1.4,0,0) (open symbols) as a 
function of magnetic field measured on IN14 (circles) and FLEX (squares). The solid line is a fit 
to Eq. (|5|) and the dotted line is a linear fit. 
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